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Abstract 

We study the effects of quantum fluctuations on excitation spectra in the two-dimensional Heisen- 
berg antiferromagnet by means of the 1/S expansion. We calculate the spin- wave dispersion and 
the transverse dynamical structure factor up to the second order of 1/5 in comparison with in- 
elastic neutron scattering experiments. The spin-wave energy at momentum (tt, 0) is found to be 
about 2% smaller than that at (jr/2,n/2) due to the second-order correction. In addition, we study 
the dimensional crossover from two dimensions to one dimension by weakening exchange couplings 
in one direction. It is found that the second-order correction becomes large with approaching the 
quasi-one dimensional situation and makes the spin-wave energy approach to the des Cloizeaux- 
Pearson boundary for S = 1/2. The transverse dynamical structure factor is also calculated up to 
the second order of 1/S. It is shown that the intensity of spin- wave peak is strongly reduced while 
the intensity of three-spin-wave continuum becomes large and exceeds that of the spin-wave peak 
in the quasi-one dimensional situation. 

PACS numbers: 75.40. Gb, 75.10.Jm, 75.30.Ds 
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I. INTRODUCTION 



The physics of low- dimensional spin systems has attracted much interest for past decades. 
Although the quantum fluctuation is expected to be large in two dimensions, there ex- 
ist strong evidences that the quantum Heisenberg antiferromagnet (QHAF) with nearest- 
neighbor coupling in a square lattice exhibits the Neel long-range order at zero temperature.— 
Under the presence of the long-range order, the linear spin-wave (LSW) theory works rather 
well, 2 ' 3 A natural way to include quantum fluctuations is an expansion in terms of 1/S, 
where S is the magnitude of spin, because the LSW theory is made up of a leading-order in 
the 1/S expansion. Such attempts have been done and turned out to be useful . 4 ' 5 ^ 7 ' 8 ' 9 

In our previous paper,— basing on the Holstein-Primakoff transformation,— we calcu- 
lated corrections up to the second order of 1/S* in various physical quantities such as the 
spin-wave dispersion, the sublattice magnetization, the perpendicular susceptibility, and the 
spin-stiffness constant. We also calculated the dynamical structure factors of the transverse 
and the longitudinal components up to the second order of 1/S.— At that time, however, 
available experimental data of inelastic neutron scattering (INS) were limited to the momen- 
tum transfer in a narrow region of the Brillouin zone (BZ). Therefore, our study was just 
a demonstration of usefulness of the 1/S expansion. Now that INS experiments provide us 
the information of the excitation spectra in the whole BZ, it may be interesting to calculate 
the excitation spectra in the whole BZ in comparison with recent experiments. We show 
that the second-order correction makes the spin wave energy at momentum (jr, 0) about 2% 
smaller than that at (tt/2,tt/2). This difference is smaller than the value 7 ~ 9% obtained 
by the series expansio n 11 ! 12 and the Monte Carlo simulation.— We have so far not been 
able to find why the 1/S* expansion within the second order gives different results, since the 
higher-order corrections is expected to be quite small. Note that the spin-wave dispersion 
was recently measured by the INS experiment for Cu(DCOO)2-4D 2 (CFTD), revealing the 
6% difference . 14 ^ This material is believed to be well described by the S = 1/2 Heisenberg 
model within the nearest-neighbor coupling.— In addition to the spin-wave energy, we cal- 
culate the transverse dynamical structure factor up to the second order of 1/S. It consists 
of the 5-function-like peak of one spin-wave excitation and the continuum of three spin-wave 
excitations. The second-order correction is found quite small in the spin- wave-peak inten- 
sity due to a cancellation of various second-order processes. The result is compared with 
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the recent experiment^ as well as other l/S'-expansion study 17 based on the Dyson-Maleev 
transformation and the series expansion studies . 11 ' 12 In the three-spin-wave continuum, such 
a cancellation in the second-order processes is not severe, and the substantial intensities 
come out. This is consistent with our previous study»& and others . 11 ' 12 ' 17 

Another purpose of this paper is to study the crossover behavior in the QHAF from two 
dimensions to one dimension with weakening exchange coupling in one direction. In purely 
one dimension, of course, the antiferromagnetic long-range order disappears due to quantum 
fluctuation and therefore the concept of spin waves breaks down. Carrying out the 1/S 
expansion, we demonstrate that the second-order corrections increase with approaching the 
quasi-one dimensional situation. The second-order correction works to increase considerably 
the sublattice magnetization, although the first-order correction makes it decrease, in the 
quasi-one dimensional situation. Interestingly, the spin-wave dispersion is found to approach 
to the curve known as the des Cloizaux- Pearson boundary in the S = 1/2 QHAF chain.— 
At the same time, for the spin-wave peak in the transverse dynamical structure factor, the 
first-order correction makes the intensity decrease but the second-order correction makes it 
increase. In the quasi-one dimensional situation, the former is much larger than the latter, 
and the net intensity is strongly reduced from the LSW value. On the other hand, the 
intensity of three-spin-wave continuum by the second-order correction increases and exceeds 
the spin-wave-peak intensities in the quasi-one dimensional situation. This contrasts with the 
description of using spinon^ to describe the large intensity of the spectral continuum. The 
above three characteristics in the quasi-one dimension, (a) the spin-wave energy approaches 
to the des Cloizaux-Pearson boundary, (b) the spin-wave-peak intensity decreases, and (c) 
the intensity of three-spin-wave continuum increases, are suggesting a close relation to the 
purely-one dimensional behavior that the spectra are described by continuum of two spinons. 

The present paper is organized as follows. In Sec. m the Heisenberg Hamiltonian is 
expressed in terms of the 1/S-expansion. The Green's functions for spin waves are introduced 
in Sec. IIHI The sublattice magnetization is calculated with the help of the Green's functions 
in Sec. IIVI The spin-wave dispersion is calculated in Sec. EJ and the transverse dynamical 
structure factor is calculated in Sec. IVII Section rvTTl is devoted to the concluding remarks. 
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II. HAMILTON! AN 



We consider the Heisenberg Hamiltonian on the square lattice with directional anisotropy 
of exchange couplings: 

H = J ■ Se+ a + J' ^ Sg ■ S^ + b, (2.1) 

t t 

where i runs over all lattice sites and I + a and i + b indicate the nearest neighbors to the 
£th site in the positive x and y directions, respectively. Quasi-one dimensional situations 
are realized by weakening the exchange coupling J' in the y direction. 

Introducing the Holstein-Primakoff transformation,— we express the spin operators in 
terms of boson annihilation operators dj and bj (and their Hermite conjugates), 

S* = S-a\a u (2.2) 

St = (Sr)t = v^S/iOSH (2-3) 

S* = -S + b]bj, (2.4) 

S+ = (Sr)t = V2Sb]f 3 (S), (2.5) 

where 

with = a\di and Indices i and j refer to sites on the "up" and "down" sublattices, 

respectively. Substituting Eqs. ()2.2|) - ()2.5|) into Eq. (j2.1j) we expand the Hamiltonian in 

powers of 1/S as 

H = -S 2 N(J + f) +H + H 1 + H 2 + ••• , (2.7) 

with N the number of lattice sites. 
The leading term Hq is expressed as 

H = JS y^qtqj + 2b i+a b i+a + aA+a + aA-a + 44+* + Q ? b I-a) 

i 

+ J'S y^(2gtqj + 2b i+b b i+h + aA+b + aA-b + + a ! fo !-b)- ( 2 - 8 ) 

i 

We diagonalize Hq by rewriting the boson operators in the momentum space as 

/ 2 \ 1/2 

a * = ( ]v ) 5Z flkexp ^ k ' ( 2 - 9 ) 
/ 2 \ 1/2 

6 i = (ivj E 6kex P^ k - r ^' ( 2 - 10 ) 



and by introducing the Bogoliubov transformation, 



<4 = 4«k + m k /3_ k , 6_ k = m k al + 4/3-k, 



where 



with 



l + e k 
L 2e u 



1/2 



m k = -sgn(7 k ) 



2e 



k J 



1/2 



EE -X k 4 



Ck 

7k 

c 



I 1 - 7k) , 

(cosA; x + (cosk y )/(l + C), 
J' I J. 



(2.11) 
(2.12) 

(2.13) 
(2.14) 
(2.15) 



Momentum k is defined in the first magnetic Brillouin zone (BZ). The sgn(7 k ) denotes the 
sign of 7 k , which is absorbed into the definition of x k . For the study of the isotropic exchange 
coupling (( = 1), we have neglected this factor because 7 k is always positive in the first BZ. 
For the anisotropic coupling, however, this redefinition of x k is necessary, because 7 k is 
negative in a certain region of the first BZ. After this transformation, we have 



tf = 2JS(1 + C)5>k-1) 

k 

+ 2J5(l + C)^ekK«k + /5 k /3 k ) 



(2.16) 



This expression is the same as that for the isotropic coupling, except for the first factor 
2J5(l + 0- 

The first-order term Hi can be expressed in terms of spin-wave operators through the 
same procedure as above. The result for the anisotropic coupling is given by the previous 
expression in Ref. with simply replacing JSz by 2JS(1 + 77): 



+ 2J l ( !+° ^^(1 + 2-3- «l 4 

1234 



X 



2SN 



+ (2a{(3„ 2 a 3 a 4 B$ u + 2[3l,(3_ 1 p_ 2 a 3 B{%, + a\a\^_ 3 ^BfL + H.c.) 



(2.17) 



5 



with 

k 

Momenta ki, k2, k3, • ■ ■ are abbreviated as 1, 2, 3, • ■ ■ . The first term arises from setting 
the products of four boson operators into normal product forms with respect to spin-wave 
operators. The second term in Eq. (j2.17j) represents the scattering of spin waves. The 
Kronecker delta <5g(1 + 2 — 3 — 4) represents the conservation of momenta within a reciprocal 
lattice vector G. The vertex functions B«'s in a symmetric parameterization are the same 
as those given by Eqs. (2.16)-(2.20) in Ref. y, so that they are omitted here. 

The second-order term H 2 is composed of products of six boson operators. Writing it in 
a normal product form with respect to spin-wave operators, we have 

H 2 = 2 J ^ } t ~ E [^(^(4^ + + C 2 (k)(4/?l k + /3_ k a k ) + ■■•]. (2.19) 
Neglected terms are unnecessary for calculating corrections up to the second order. The 
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Note that 



explicit forms of Ci(k) and C^k), are given by Eqs. (2.22) and (2.23) in Ref. 
Ci(k) and C^k) diverge as l/e k with |k| — > 0. 

III. GREEN'S FUNCTION 

For systematically carrying out the 1/S-expansion, it is convenient to introduce the 
Green's functions for spin-waves, 

G aa (k,t) = -i(r(a k (t)4(0))>, (3.1) 

G af3 (k,t) = -z(T(a k (t)/l_ k (0))>, (3.2) 

G Pa (k,t) = -i{T(pi k (t)al(0))), (3.3) 

G pp Qt,t) = -*(T(/?l k (t)/3_ k (0))>, (3.4) 

where (• • • ) denotes the average over the ground state, and T is the time-ordering operator. 

In this paper, we measure energies in units of 2JS(1 + £). The unperturbed propagators 
corresponding to H are given by 

Gl(k,c) = [u-e^ + iS]- 1 , (3.5) 
G%(k,u) = G° a (k,^) = 0, (3.6) 
G%(k,u) = [-cu-e k + t5]- 1 . (3.7) 
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Sk iP ,q,[k h p q] S kpq[ktp q] fik,p,q,[k*p-q] 

(b) 




FIG. 1: Second-order diagrams for the self-energy, (a) £ Qa (k, u) and (b) E aj g(k, uj). Solid lines 
represent the unperturbed Green's functions. 



The self-energy is defined by 

G^(k,uj) = G%(k,iu) + £G^,(k,w)£^(k,w)GW(k, 

ix'v' 

It is expanded in powers of 1/(2S), 

£^(k,u,) = ^E«(k,w) + -ipSgCk,^) + • • • . 
^From Hi we have the first-order correction as 

E<2(k,o;) = S«( k '^) = ^, E«(k,w) = Eg(k, W ) = 0. 



(3.8) 



(3.9) 



(3.10) 



( 2} 

The second-order term S^ (k, uj) is obtained from the second-order perturbation, whose 
diagrams are shown in Fig. ^ We obtain formally the same expression for the self-energy 



as in our previous paper: 
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sL 2 2(k,^) 



Ci(k) 



E 2 ^ 



B 



2 p2 p2 p2 

pq 

(4) 

k,p,q,[k+p-q] 



B 



(6) 

k,p,q,[k+p-q] 



uj — e 



p c q 



ek+p-q + iS uj + e p + e q + e k+p _ q - z<5 



(3.11) 



Sg(k,^) = Sg(-k,-c) 

= C 2 (k)+ (^EW +p _ q s g n( 7G ) 



x B 



(i) 



,5. 



(6) 



2(e p + e q + e k+p _ c 



k,p,q,[k+p-q] k,p,q,[k+p-q] ^2 _ _|_ 6q _|_ e k+p _ q ) 2 + Z<$ ' 



(3.12) 
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where 5 — > 0, and [k + p — q] is the vector k + p — q reduced to the 1st BZ by a reciprocal 
vector G. We have used the relations 

S [k+p-q],q,p,k = S g n (^G)^ iPiqi[k+ p_ q ], ^ 
5 q 6 [k+p-q],k,p = S g n (7G)5^ q[k+p _ q] . 

The terms divergent with k — > in Ci(k) and C^k) are canceled by the second-order 
perturbation terms in Eqs. (|3.11|) and (|3.12j) . One can prove £/J(k — > 0,u; = 0) — ► from 
these equations. 



IV. SUBLATTICE MAGNETIZATION 



Once the Green's function is known, the sublattice magnetization is calculated from the 
relation, 



M = S— (a\ ai ) 

2 E, lim + r°°^ie^Ul G ^M 



S 



N 



2tt 



with r\ — > + . After carrying out the integration with respect to u, we obtains 

M 2 



Af = 5 - AS 



(2sy 



with 



AS 



Mo 



2 

N 



l E ^E2(k,- £k ) 



■ J 4^jc , 



2 



^ ' pq 



?2+ m 2 ) | jB (6) 



k,p,q,[k+p-q] 



+ 



24m k sgn(7 G )5 ; 



(ek + e p + e q + e k+p _ q ) 

>(6) 



(4) B V 
k,p,q,[k+p q] k,p,q,[k+p-q] 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



4 - ( e P + e q + e k+p _ q ) 2 

Here [k + p — q] was defined before, and G = k + p — q — [k + p — q] . The zeroth-order 
correction AS represents the well-known "zero-point" reduction in the LSW theory. To 
evaluate Eq. (|4.4|) . we sum up the values of points of k in the 1/4 part of the first BZ 



S 



TABLE I: Sublattice magnetization 



J'/J 


1 


0.5 


0.1 


0.075 


0.05 


AS 


0.196 


0.213 


0.355 


0.391 


0.445 


M 2 


0.0035 


0.024 


0.323 


0.481 


0.818 



and Nl points of p and q in the first BZ, with Nl = 20, 48. For J' / J = 1, the convergence 
is very good; M 2 = 0.0035059 for N L = 20, and M 2 = 0.0035065 for Nl = 48. 

Table|l]lists the values of AS and M 2 for several values of J' / J. These values are evaluated 
for Nl = 48. For the isotropic coupling, we reproduce the values obtained previously.-^ The 
zero-point reduction AS* increases with decreasing values of J'/J. On the other hand, M 2 is 
found always positive, tending to cancel the zero-point reduction. The value increases with 
decreasing values of J'/J, and becomes comparable to AS around J'/J = 0.1, suggesting 
the applicability limit of the expansion for the case of S = 1/2. 



V. SPIN- WAVE DISPERSION 



Within the second order in 1/5", the renormalized spin-wave energy e k in units of 2JS(1 + 
() is obtained from 



e k = e k + — Ae k + 7 — ^(k, e k ). (5.1) 



— A 1 

2S k+ (2Sy 

From this equation, we define the renormalized spin-wave velocity V x along the x direction by 
V x = lim^^o 2 JS(1 + C)<W (fikx) with k y = 0. Thus the renormalization factor is expressed 
as 

v ~ 2JS(1 + C) 1/2 ~ 2S + (2S)^ [ ] 



with v i = A given by Eq. (|2.18|) . In the following numerical evaluation, we divide the first 
BZ into N? meshes with Nl = 64. 



A. Isotropic case 

Figure 121 shows the spin-wave energy 2JS(1 + Qe^ as a function of momentum for the 
isotropic coupling {J'/J = 1) with S = 1/2, in comparison with the experimental data 
taken from the INS for Cu(DCOO) 2 -4D 2 (CFTD) . 14 i 15 Momentum is measured in units of 
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(71,0) (7[/2,7l/2) 



FIG. 2: Spin-wave energy as a function of momentum for the isotropic coupling (J' / J = 1). 
S = 1/2. The dotted, broken and solid lines represent the values calculated within the LSW theory, 
up to the first-order correction, and up to the second-order correction, respectively. Experimental 
data are taken from Ref. Ila . Inset indicates high symmetry lines which momentum varies along. 



(nearest neighbor distance) -1 . In the whole BZ, both the first and second order corrections 
make the spin-wave energy larger. The curve along the line (0, 0) — (jr, 0) has already been 
reported in our previous paper.— The dispersion along (vr/2, tt/2) — (n, 0) is completely flat 
within the first-order correction. The second-order correction makes the excitation energy 
at (7r, 0) about 2% smaller than the energy at (7r/2,7r/2). Explicitly they are e( 7r /2, 7r /2) = 
1.196 and e^o) = 1-179. A previous series expansion study predicted the energy difference 
about 7%,— and a recent study gave about 9% difference, that is, = 1-192 and 

e^o) = 1.09.— A Monte Carlo simulation has given e( 7r /2, 7r /2) = 1-195 and £( n ,o) = 1-08. 13 
These values at (7r/2,7r/2) agree well with our value, while the values at (tt, 0) is rather 
different from our estimate. The experimental data indicate that the excitation energy at 
(vr,0) is 6% smaller than that at (vr/2, vr/2). 



B. Anisotropic case 

Figure El shows the renormalized spin-wave energy as a function of momentum along 
(0,0) — (vr,0) for J' /J < 1. As the same as the isotropic coupling, both the first-order 
and the second-order corrections are found to be positive, making the energy larger. Both 
corrections increase with decreasing values of J'/ J. In quite small interchain couplings 
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1 



1 

I of 



» * 1 

'to 




(0,0) (71,0) 

FIG. 3: Spin-wave energy as a function of momentum for anisotropic couplings, (a) £ = 0.5, (b) 
£ = 0.1, (c) C = 0.075, and (d) £ = 0.05. The dotted, broken and solid lines represent the values 
calculated within the LSW theory, up to the first-order correction, and up to the second-order 
correction, respectively. The thin solid line labelled "ID Exact" represents the des Cloizeaux- 
Pearson boundary. 



TABLE II: Renormalization of spin-wave velocity 



J'/J 


1 


0.5 


0.1 


0.075 


0.05 




0.158 


0.174 


0.272 


0.287 


0.306 


v 2 


0.021 


0.053 


0.130 


0.141 


0.155 



(J'/J < 0.1), the excitation energy seems approaching the des Cloizeaux-Pearson boundary 
in one dimension.— As shown in Table |HJ the renormalization constant Z v seems approaching 
7r/2, corresponding to the value of the boundary. 




11 



VI. DYNAMICAL STRUCTURE FACTOR 



The dynamical structure factor is an important quantity, since it is directly related to 
the INS spectra. We have already reported the expression within the second order of 1/5* in 
the isotropic coupling situation.— As is evident from the forms of H and Hi, the formulas 
for the anisotropic coupling are formally the same as those for the isotropic coupling. 

We consider only the transverse component, which is defined by 

S^{k,u) = ^J dte^(Q u(s) (M)Q u(s) (k,0)t), (6.1) 



where 



with 



Q n(s) (k) = S+(k)±S+(k), (6.2) 



/ 2 \ 1/z 

S a + (k) = K(k)]t = (-) £ S+expHk-r,), 

' i 

S*(k) = [5 6 "(k)]t= (ly^S+expi-ik-rj). 



(6.3) 
(6.4) 



We need the "uniform" and "staggered" parts because the momentum is defined inside the 
first BZ. They are labelled as the suffix "u" and "s", and correspond to upper and lower 
signs in Eq. (|6.2|) . respectively. 

We start by introducing the operators, 

Y a + (V = [17(k)]t = [4S a +(k) - ro k #(k)]/(2S) 1/a , (6.5) 
F+(k) = [iy(k)]t = [-ro k 5+(k) + 45+(k)]/(2S) 1 / 2 , (6.6) 

and the associated Green's functions, 

/oo 
dte^(T[Y^(k,t)Y-(k,0)]). (6.7) 
-oo 

Then, with the help of the fluctuation-dissipation theorem, we have 



5+- ) (k,^) = 25(4±m k ) 2 (~ 



xIm[F aa (k, u) ± F a/3 (k, u) ± F^ a (k, u) + Fp p (k, u)), 

where the upper (lower) signs correspond to the uniform (staggered) part. 
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(6.8) 



For calculating F^ u (k, u), we expand the operator Y+(k) in terms of spin- wave operators 
with the help of the HP transformation and the Bogoliubov transformation. After lengthy 
calculations, we have 

234 

x (M$ u (3_ 2 a 3 a 4 + M^ u 4^_ 4 + •••), (6.9) 



Y+(k) = D(3\ 
1 2 



^5 G (k + 2 - 3 - 4)^W^4Sgn(7 G ) 



2S 1 iV 

234 

x (M( 2 2 ) 34 /3_ 2 a 3 a4 + M^a^/?^ + ••■), (6.10) 



where 



fl = 1 _Ag_lAg(l+jAg), (6 . n) 



2S 4 (25) 2 

'k234 " ~~ ' 

f (2) 
'k234 

with G = k + 2 — 3 — 4. The first-order and second-order corrections in Eq. ()6.11j) arises from 



M k234 = ~^2 + sgn(7 G )x k x 3 X4, (6.12) 

(2) 

M k234 = X 3^4 - sgn(7 G )x k x 2 . (6.13) 



setting four and six boson operators in the HP transformation into the normal product forms 
with spin-wave operators, respectively. Thereby, the second terms in Eqs. ()6.9|) and ()6.10|) 
are normally ordered. Note that sgn(7 G ) arises from the phase difference in the definitions, 
Eqs. (JEHSJ) and (Q- 

With the use of Eqs. ()6.9j) and ()6.1()j) . F^ u (k,u) is expanded up to the second order by 
the diagrams shown in Fig. |3J Explicitly, it is given by 

F, u (k,u) = D 2 Gl u (k,u)5^ + Gl(k,u)^ 2 ^XKuj)Gl(k,uj) 

+ I^(k,u;)Gl(k,u) + G^(k,^(k,^) + J^(k,u), (6.14) 

where each term in Eq. ()6.14j) corresponds to the diagrams (a), (b), (c), (d), and (e), respec- 
tively. Explicit expressions for /^(k, u), I^ik, u), and J^ u (k,uj) are given by Eqs. (4.5)- 
(4.22) in Ref. H 

The dynamical structure factor is obtained by substituting Eq. ()6.14|) into Eq. (j6.8J) . It 
consists of the 5-function-like peak of the one spin-wave excitation and the continuum of 
three spin-wave excitations: 

Sfofr") = P { uU k ^ - ek ) + ( 6 - 15 ) 
13 



FIG. 4: Diagrams for .F Mi ,(k, uj). Solid lines represent the unperturbed Green's functions Gr* (k, u). 



For the term of one spin-wave excitation, the bare energy e k is to be replaced by the renor- 
malized value e k given by Eq. (J5.1j) . However, the spectral weight p^(k) within the second 
order of 1/5 is safely evaluated by putting to = e k in Eq. ()6.14j) . It is expressed as 



with 



pi 1 ( ) s) (k)=25(4±m k ) 2 



d 



u{s),l 

~2S~ 



+ 



d 



u(s),2 



d u (s),i 

d u (s),2 



-2AS. 
1 

~2 



AS(l + A£-)=F-£$(k,e k ) 



4) 5> 



2 p2 p2 p2 
k^p^q^k+p-q 



pq 



-\B {4) I 2 

1 k,p,q,[k+p-q]l 



+ 



l B (S) 12 
1 k,p,q,[k+p-q]l 



2 



Ck+p-q) 2 (fk + e P + Cq + Ck+p-q) 2 

+ (1) ' £ 2 ^<Cp-q(<L, [k+P - q] ± sgn( 7G )ii41 q , [k+p _^ 

sgn( 7G ) J B k 6 ^ q ^ 



pq 

B (4) 

k,p,q,[k+p-q] 



e k+p-q; 



p qj 



ek 



^k+p-q 



(6.16) 



(6.17) 



(6.18) 



The upper (lower) signs correspond to the uniform (staggered) part. The first-order correc- 
tion Eq. ()6.17|) arises from the first term of Eq. ()6.14|) . In the second-order correction given 
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by Eq. ([6.18)1 . the first term arises from the first term of Eq. (|6.14|) . and the second term 
arises from the second term of Eq. ()6.14j) , 



(2S')2^ a/3Vrv ' uy; ^ /3/3 
The third term of Eq. ()6.18|) is equivalent to 

1 <9£i 2 2(k,cu 



(2S) 



(2S) 2 Ooj 

and arises from the second term of Eq. ()6.14j) , 



G° aa (k,u) -j- Eg (k, w ) (k, u,) . 



This is related to the second-order correction to the residue of the spin-wave pole in 
G Qa (k,cj), 

1 1 dZ ( al{k,LJ 



1 1 


SE^(k,aO 




1 (25) 2 


eta 


w=e k 



1 + 



(25) 2 diu 



The fourth term arises from the third and fourth terms of Eq. ()6.14|) , 

V(k, cu)G° uu (k, u) + G°(k, u)I^(k, u). 



No contribution arises from the last term of Eq. ()6.14j) . J MJy (k, tu). Note that the main 
momentum dependence around k = (0, 0) arises from the prefactors, (4 + ^k) 2 oc e k and 
(4 - "ik) 2 oc l/e k . 

The three-spin-wave continuum arises only from the second-order corrections. Note that 
the first term of Eq. ()6.14|) has no contribution. After careful evaluation of other terms in 
Eq. ()6.14|) . we obtain 

2 



p%(k,u) = 2S(4±m k ) 2 ^(£) 



to — e r 



e k+p-q; 



A 2 k p °i k +p _c i 
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k,p,q,[k+p-q] 



2sgn( 7G )i? k 6 p , qJk+p 

Ck + Cp + Cq + e k+p-q 



£k - e t 



p q 



k,p,q,[k+p-q] 
1 2 



(6.19) 



The spectral shape may be modified by the renormalization of spin-wave energies and by 
taking account of scattering spin waves due to mutual interaction, which terms are present 
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in Hi. Therefore, it may be difficult to determine the spectral shape in a consistent way 
with the 1/5* expansion. However, the total intensity, which is given by 

POD 

O k )= / d ^i(l)( k ^)> ( 6 - 2 °) 
may be safely evaluated from Eq. ([6. 19)1 . Note that around k = (0, 0) 



S, q ,[k + p- q] ± sgn(7G)M( 2 ) q [k+p _ q] cx e k , (6.21) 



e k — e p — e q — e k+p-q e k + e p + e q + e k+p-q 

?( 4 ) -r am .f„ \n( 6 ) 



B k,p,q,[k+p-q] =F S g n (7G)-B kp)q [k+p _ q] 



oc e k , (6.22) 

e p + e q + e k+p-q 

and (ik+rrik) 2 ^ oc const., (4 _ Wk) 2 4 ^ l/ e k; we n °ti ce the dependences around k = (0, 0) 
as Iu (k) oc e k , li (k) oc const. 

In the numerical evaluation of Eqs. ()6.18|) and ([6. 19)1 . we sum up the values on iV| points 
of p and q in the first BZ, with Nl = 64. 



A. Isotropic case 

Figure El shows the spin- wave-peak intensity and the intensity of three-spin- wave contin- 
uum as a function of momentum along high-symmetry lines for S = 1/2. Using the extended 
zone scheme, we assign the staggered part for line (0,0) — (7r/2, 7r/2) to the values for line 
(7r,7r) — (7r/2,7r/2) and also the staggered part for line (0,0) — (jr, 0) to the vales for line 
(7r, 7r) — (tt, 0). The uniform part is assigned inside the first BZ. At the zone boundary of the 
reduced BZ, the uniform and staggered parts coincide with each other. The second-order 
corrections to the spin-wave-peak intensity becomes one order of magnitude smaller than the 
first-order correction, due to a cancellation among contributions of four terms in Eq. (|6.18j) . 
The intensity is almost determined within the first-order correction. Thus the correction rel- 
ative to the zero-th order value is independent of momentum. We obtain around k = (0, 0), 

p«(k)=0.215|k|, p«(k) = 1.72/|k| (S = l/2). (6.23) 

These values should be compared with the l/S'-expansion analysis based on the Dyson- 
Maleev transformation,— 0.202 and 1.86. A series expansion analysis by Singh^S gives the 
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FIG. 5: Transverse dynamical structure factor as a function of momentum for the isotropic cou- 
pling. S = 1/2. The dotted, broken and solid lines represent the intensity of one spin- wave excita- 
tion calculated within the LSW theory, the first-order correction, and the second-order correction, 
respectively. Thin solid line represents the intensity of three spin-wave excitations. Experimental 
data are taken from Ref. Ila . Inset indicates high symmetry lines which momentum varies along. 



values, 0.246 and 2.10, while a recent analysis by Zheng et alM gives the values, 0.216 and 
1.86. 

The small second-order correction to p^K (k) does not necessarily mean small three-spin- 
wave continuum. At the zone boundary (vr, 0), for example, we have l^(ir,0) = 0.143 in 
addition to p^(tt, 0) = 0.618. Such a considerable intensity of three spin-wave continuum 
has been predicted in our previous paper— and others . 12 > 17 It varies as proportional to e£ at 
the zone center, and converges to a constant value at (7r, 7r). 



B. Anisotropic case 

Figure |H1 shows the spin- wave-peak intensity and the intensity of three-spin- wave contin- 
uum along the symmetry line (0, 0) — (tt, 0) for several anisotropic couplings. 5* = 1/2. 

The spin-wave-peak intensity is reduced by the first-order correction, but is increased by 
the second-order correction. Note that, although the residue of the spin- wave pole in the 
Green's function is reduced by the self-energy, the other terms of the second-order correction 
work to increase the intensity. The second-order correction increases with decreasing values 
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FIG. 6: Transverse dynamical structure factor as a function of momentum along the symmetric 
line (0,0) - (vr,0) for anisotropic couplings; (a) J'/ J = 0.5, (b)J'/J = 0.1, (c)J'/J = 0.075, and 
(d) J' / J = 0.05. S = 1/2. The dotted, broken and solid lines represent the spin-wave-peak intensity 
calculated within the LSW theory, up to the first-order correction, and up to the second-order 
correction, respectively. The thin solid line represents the three-spin-wave continuum intensity. 

of J'/ J, but the reduction due to the first-order correction is much larger than the gain due 
to the second-order correction. As a result, the spin- wave-peak intensity is strongly reduced. 

On the other hand, the intensity of three-spin-wave continuum increases with decreasing 
values of J' /J. It exceeds the spin- wave-peak intensity in the quasi-one dimensional sit- 
uation. Such large intensities of continuum spectra have been observed in the recent INS 
experiments on the quasi-one dimensional QHAF such as KCUF3 21 and BaCu2Si207.— 

VII. CONCLUDING REMARKS 

We have systematically carried out the 1/S expansion up to the second order on the basis 
of the HP transformation in the two-dimensional QHAF. We have calculated the spin-wave 
energy in the whole BZ, in comparison with the recent INS experiment for CFTD . 14 ^ We 
have found that the spin-wave energy at (vr,0) is about 2% smaller than that at (7r/2,7r/2) 
due to the second-order correction. This is a correct tendency, but the value is somewhat 
smaller than the experimental value 6% and other theoretical estimates 7 — 9% . 1:L i 12 i 13 We 



18 



have so far not been able to find the reason for the difference, since the corrections higher 
than the second order of 1/S is expected to be quite small. We have also calculated the 
transverse dynamical structure factor. The second-order correction is found extremely small 
in the one-spin-wave-peak intensity due to the cancellation in the second-order terms, while 
it gives rise to substantial intensities of three-spin-wave continuum. This is consistent with 
previous studies.— i n i 12 i 17 

Canali and Wallin^ reported that the value of the perpendicular susceptibility x± is 
different from the value in our previous paper 5 (for J' /J = 1), although they confirmed 
the same values for the second-order correction to the sublattice magnetization M2 and 
for the spin- wave velocity V x . The difference is not due to numerical errors, since the 
convergence with respect to Nl has been carefully checked with changing Nl = 160, 320, 480 
(see Ref. 0). As already discussed there, the value in Ref. satisfies the hydrodynamic 
relation, V x = (ps/x) 1 ^ 2 i n an appropriate unit, with independently-evaluated spin-stiffness 
constant p s . As regards the dynamical structure factor, the difference between the values 
by Canali and Wallin and the present values might have the same origin as the difference in 
the perpendicular susceptibility, because the diagrams for the perpendicular susceptibility 
are closely related to those for the transverse dynamical structure factor. 

The second-order correction is expected to become more important in quasi-one dimen- 
sional systems. We have studied the crossover from two dimensions to one dimension by 
weakening the exchange coupling in one direction. All formulas are found formally the same 
as those for the isotropic coupling with replacing JSz by 2JS(1 + (). It is found that the 
excitation energy is pushed up by the first-order and second-order corrections. With ap- 
proaching the quasi-one dimensional situation, the corrections make the excitation energy 
close to the des Cloizeaux-Pearson boundary in the one-dimensional QHAF for S = 1/2. 
In the transverse dynamical structure factor, the spin- wave-peak intensity is reduced by the 
first-order correction, but is increased by the second-order correction. The former exceeds 
the latter in the quasi-one dimensional situation, and thereby the peak intensity is strongly 
reduced from the LSW value. On the other hand, the intensity of three-spin-wave continuum 
is found to become larger and exceeds the spin-wave-peak intensity. 

In purely one-dimension, spin-one excitations are considered excitations of two spinons 
of a spin-one-half excitation. In this respect, the spin wave might be considered as a bound 
state of two spinons. Our finding that the weight pi(k) of the spin- wave peak decreases with 
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J'/ J — > is consistent with this picture. Large intensities of three-spin-wave continuum 
might be replaced by two-spinon continuum in the one- dimensional limit. Recently, INS 
experiments have been carried out at low temperatures in the quasi-one dimensional systems 
such as KCuF 3 21 and BaCu2Si207,— and large broad spectra have been observed in addition 
to a peak in the transverse dynamical structure factors. This behavior as well as the behavior 
of the longitudinal component have been analyzed by the chain-mean-field and random phase 
approximation . 23 ! 24 It may be interesting to analyze these data in terms of the 1/S expansion 
by starting from a detailed three-dimensional model with directional anisotropy. 
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